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The development of the microbunching instability is studied for the linac of the proposed Shanghai 
Soft X-ray Free-Electron-Laser facility (SXFEL) by analytic formulae as well as numerical simula- 
tions with the aid of two different codes. The process is investigated in detail and the growth rates 
(gains) of the instability under various conditions are compared. The results indicate that the limi- 
tations from numerical computations in the present simulation model must be taken into account. 
Moreover, the appearance of higher-order mode excitations in the simulations suggests that further 
improvement of the current theory is necessary. A mechanism of introducing shot noise into beam 
profile as the beam passes through a chicane is proposed. The key issues that drive the instability 
are analyzed. 
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I. INTRODUCTION 

The microbunching instability in the linac of a free- 
electron Laser (FEL) facility has always been a problem 
that degrades the quality of electron beams. The insta- 
bility is driven by various effects, such as the longitu- 
dinal space charge (LSC) [l(, coherent synchrotron ra- 
diation (CSR) 0,3) and linac wakefields. As the beam 
passes through a bunch compressor, (e.g., a magnetic chi- 
cane) the energy modulation introduced by those effects 
is transformed into density modulation and thus the in- 
stability develops. In a FEL facility, there is usually more 
than one bunch compressor and the overall growth of the 
instability is the product of all the gains in each com- 
pressor. As a result, the final gain of the instability can 
become significant. On the other hand, the FEL process 
has a high demand for electron beam quality in terms of 
peak current, emittance, energy spread, etc. Therefore 
without effective control, the microbunching instability 
can damage the beam quality so seriously that the whole 
FEL facility fails. Thus learning how to control and/or 
reduce the instability is a key to the success of a FEL 
project. 

The Shanghai Soft X-ray FEL Facility (SXFEL) 
project has been proposed and the feasibility study has 
been finished. The facility is expected to be built in a few 
years. Once it is built, it will be the first X-ray FEL fa- 
cility in China. It is a cascading high-gain harmonic gen- 
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eration (HGHG) FEL facility working at the wavelength 
of 9 nm in the soft- X-ray region. The basic beam pa- 
rameters at the exit of linac are: beam energy 840 MeV, 
peak current 600 A, and normalized transverse emittance 
less than 2.0 mm.mrad. Both analytical and numerical 
studies show that the microbunching-instability-induced 
growth of the global and slice energy spread in the linac is 
not ignorable without proper control. Moreover, the in- 
stability reduces the smoothness of the longitudinal beam 
current profile as well. One way to control the instability 
is to increase the uncorrelated energy spread of the beam 
by a laser heater [J], which will be implemented in the 
SXFEL. 

In this article, we study the microbunching instability 
based on the design parameters of the SXFEL in order 
to gain better understanding of the instability process so 
as to uncover the limitations in analytic formulae and 
numerical simulations. We realize that the higher har- 
monics of beam current can be excited in a bunch com- 
pressor which may lead to modification to the linear the- 
ory of the instability gain. Numerical parameters, such 
as the half-width and the order of the Savitzky-Golay 
filter [5[ employed in the computation of the longitudi- 
nal space-charge (LSC) force and coherent-synchrotron- 
radiation (CSR) effect, play an important role in the 
energy spread of the local slice. Moreover, the shot 
noise introduced inevitably in a bunch compressor due 
to transverse-longitudinal coupling will be important in 
a laser heater. In Sec. II, we introduce the basic struc- 
ture of the SXFEL linac, the fundamental theory of mi- 
crobunching instability, and the gain curves of the insta- 
bility in the SXFEL linac obtained by analytic formulae 
and numerical simulations based on the present design 



parameters. Comparisons of the gain curves under dif- 
ferent conditions are provided and some sensitive param- 
eters are analyzed. In Sec. Ill, some critical issues found 
in the instability study which may have significant ef- 
fects on the final results are discussed. Summaries and 
concluding remarks are given in Sec. IV. 



II. MICROBUNCHING INSTABILITY 
ANALYSIS IN SXFEL 



The basic mechanism of the microbunching instabil- 
ity driven by various effects has been studied and corre- 
sponding theories have been derived in Refs.Q], 0, |j- The 
development of the instability is similar to the amplifi- 
cation process in a klystron amplifier. The initial den- 
sity modulation and/or white noise are transformed into 
energy modulation by various kinds of impedance (e.g., 
LSC, CSR, linac wakefields) during beam transportation. 
The energy modulation accumulates and is fed-back into 
the beam after passing through a dispersive section such 
as a bunch compressor. This accumulation results in 
much stronger amplified density modulation. Addition- 
ally, the CSR effect in the dispersive section forms a pos- 
itive feedback to enhance the instability. More dispersive 
sections in a linac result in stronger instability growth. 

In the linac of the SXFEL, both S-band and C-band 
rf structures are used to accelerate the electron beam. 
One X-band structure is implemented to suppress the 
second-order nonlinear components in the longitudinal 
phase space to avoid undesired growth of transverse emit- 
tance and slice energy spread. Two magnetic chicane- 
type bunch compressors (BC1 and BC2) are used to com- 
press the beam to arrive at the required peak current. 
The layout of the SXFEL linac is shown in Fig. [TJ and 
the magnet layout and some of the Twiss parameters 
of the SXFEL linac are shown in Fig. [2] Note that S- 
band accelerating structures are used in sections LI and 
L2, whereas a C-band accelerating structure is used in 
section L3. The following discussions are based on the 
design parameters shown in Table HI which is taken from 
the SXFEL feasibility study report. Moreover, since the 
length scale in which the structural impedance is effec- 
tive is much longer than that of microbunching [6j, [7J , we 
may neglect the effects from the linac wakefields in the 
following discussions without compromising accuracy. 
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FIG. 1: (Color) Layout of the SXFEL linac. 
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FIG. 2: (Color) Magnet layout and Twiss parameters of the 
SXFEL linac. 



TABLE I: Main beam parameters used in the microbunching 
instability study for the SXFEL (obtained by ELEGANT @|, 
with laser heater turned off in simulation). 



Parameter 



bunch charge (nC) 

beam energy out of injector (MeV) 

bunch length (FWHM) at the exit of injector (ps) 

peak current before BC1 (A) 

beam radius (rms) before BC1 (mm) 

beam energy before BC1 (MeV) 

local (slice) energy spread (rms)right before BC1 (keV) 

linac length up to BC1 (m) 

R 56 of BC1 (mm) 

beam energy before BC2 (MeV) 

local (slice) energy spread (rms) right before BC2 (keV) 

linac length up to BC2 (m) 

R 56 of BC2 (mm) 

linac length after BC2 (m) 

compression ratio (BClxBC2) 

beam radius (rms) before BC2 (mm) 



Value 
~05 
130 
8 

60 
0.40 
208.4 
6.5 
17.3 
48 

422.0 
43.4 
60.3 
15 
63.7 
5x2 
0.26 



A. Microbunching instability induced by 
longitudinal space-charge (LSC) effect 

The microbunching instability for the case of linear 
compression has been discussed by Saldin et al. [l[ phe- 
nomenologically by comparing the energy distributions 
before and after compression. Consider a density mod- 
ulation at wavenumber k. Without higher harmonics of 
beam current taken into account, the gain of the instabil- 
ity driven by the wake fields upstream of the compressor 
reads [l[ 



G = Ck\R, 
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According to Ref. [0, 7 is the nominal relativistic factor 
of the electron beam with rms local energy spread c 7 in 
front of the bunch compressor, C — 1/(1 + hR^) is the 
compression ratio, h is the linear energy chirp, R^q is 
the 5-6 element of the transport matrix, Iq is the initial 



peak current of the beam, Zq = 377 f2 is the free-space 
impedance, Z tot is the overall impedance upstream of 
compressor including those of the LSC, linac wake, etc., 
and I a = 17 kA is the Alfven current. The longitudi- 
nal space-charge_impedance per unit length in free space 
takes the form 
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where rb is the radius of the beam, K\ is the first-order 
modified Bessel function of the second kind. Based on 
Eqs. ([I]) and ©, using the parameters output by the ELE- 
GANT simulation starting from a beam with ~ ±1% noise 
fluctuation in current and 1-2 keV uncorrelated (slice) 
energy spread (Fig. 01 Fig. [3]), the gains of microbunching 
instability in the region around the peak current induced 
by LSC impedance at the exits of BC1 and BC2 (first and 
second bunch compressors) are computed by the analytic 
formula and are illustrated in Fig. [3] Note that in all the 
gain curves hereafter, the gain is the pure growth "(7—1." 
The beam parameters at the exit of the SXFEL injector 
are prepared by PARMELA [fj, LLfJ] simulation employing 
one million macro-particles. The LSC impedances are 
computed separately in drift space and the accelerating 
section, 
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FIG. 3: (Color) The gains of the LSC-driven microbunching 
instability as functions of modulation wavelength at the exits 
of the first (blue) compressor BC1 and second (red) bunch 
compressor BC2 in the SXFEL linac. 

The reason why we purposely choose to use a noisy 
input instead of a smooth one is because the real beam is 
usually not ideally smooth, it always includes the ripples 
introduced by the shot noise fluctuation, the quantum 
process of field emission and the laser power jitter, etc. 
On the other hand, the linear theory is applied in the 
computation because: I. in general, the initial density 
modulation is very small, and II. the gain in the first 
compressor BC1 is small as well. Although we will point 
out in the next section that those assumptions are not 
very accurate and discrepancy exists, the linear theory 
does provide a good first estimate for the growth of the 




FIG. 4: The longitudinal current profile of the input beam. 



instability in the SXFEL linac after reviewing the results 
of numerical simulations (Fig. [7]). 
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FIG. 5: The longitudinal phase space distribution of the in- 
put beam, note that the unit of momentum is electron mass 
divided by the speed of light. 

The gain curve at the exit of BC2 appears to peak 
at higher wavelengths, because of the larger slice energy 
spread after the first compression. For this reason, we 
will see that the longer-wavelength components in beam 
current lead to stronger amplification and they dominate 
in many cases. 



B. Microbunching instability induced by coherent 

synchrotron radiation (CSR) effect in bunch 

compressor 

When a charged particle beam passes through a bunch 
compressor, coherent synchrotron radiation (CSR) can 
be emitted at wavelengths much shorter than the bunch 
itself if the density of bunch particles is modulated at 
those wavelengths. As we have already discussed, since 
the CSR effect in a bunch compressor introduces positive- 
feedback to the microbunching instability, the gain of the 
instability therefore rises rapidly. 

The gain of CSR-driven microbunching instability in 
a bunch compressor has been derived analytically [3j in 
terms of beam energy, current, emittance, energy spread 



and chirp, as well as initial lattice and chicane parame- 
ters. Assuming a beam uniform in the z-direction and 
Gaussian in transverse and in energy distributions, the 
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CSR-driven microbunching instability growth follows the 
expression 3]: 
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As described in Ref. y, the first term on right side 
of Eq. ([3]) represents the loss of microbunching in the 
limit of vanishing current, the second term (linear in cur- 
rent) provides a one-stage amplification at low current 
(low gain), and the last term (quadratic in current) cor- 
responds to the two-stage amplification at high current 
(high gain) . The functions of Fq , F\ , and F2 are defined 
in Ref. |J, and as and a x are related to, respectively, the 
local (slice) energy spread and rms transverse size of the 
beam. 

The analytical computation of the gain of the CSR- 
driven microbunching instability in the region around the 
peak current is shown in Fig. [6] In the calculation for 
the SXFEL linac, since the high-gain term in Eq. <j3j) is 
much smaller than the low-gain term (possibly due to 
the rapid increase of slice energy spread in the bunch 
compressor), we just ignore the contribution of the two- 
stage amplification. In the following section, we will show 
that this assumption is valid when comparison is made 
with simulation results. 
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FIG. 6: (Color) The gains of the CSR-driven microbunching 
instability as functions of modulation wavelength at the exits 
of the first (blue) compressor BC1 and second (red) compres- 
sor BC2 in the SXFEL linac. 



C. Discussions on final gain 

IMPACT-Z [ll| is a three-dimensional particle-in-cell 
(PIC) code designed for particle tracking in an accel- 
erator. Because it computes the electromagnetic fields 



based on fundamental equations, it may reproduce the 
dynamics of beam particles with sufficient accuracy. An- 
other reason of choosing this code is because a PIC code 
performs the smoothing of the particle distribution inher- 
ently. Moreover, it is relatively fast when a large num- 
ber of particles are tracked. The comparison of the gain 
curves at the exit of SXFEL linac obtained from ana- 
lytic formulae and from IMPACT-Z simulation is shown in 
Fig. In the IMPACT-Z simulation, for the 8-pico-second 
beam out of the injector, 512 slices are used longitudi- 
nally to compute the longitudinal space charge (LSC). 
To obtain the data points by IMPACT-Z, initial modula- 
tion of certain wavelength is added on top of the original 
beam distribution, then we track the beam through the 
linac, and the gain is computed by comparing the Fast- 
Fourier transform (FFT) patterns of the beam current 
distribution in the vicinity of the maximal current at the 
entrance and the exit of the linac. 
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FIG. 7: (Color) The gains of microbunching instability at 
the exit of the SXFEL linac section as functions of modula- 
tion wavelength computed by analytic formulae (green) and 
IMPACT-Z simulation (red) based on current machine and 
beam parameters. 

As we mentioned in the previous section, there is al- 
ways noise introduced inevitably when the electrons are 
emitted from the cathode. Therefore, in order to have a 
clearer picture about the effect of the initial noise on the 
final gain, the gain with an initially smoothed electron 
beam is computed as well. Our study shows that the un- 
correlated (slice) energy spread before the second bunch 
compressor (BC2) becomes larger for an initially noisy 



(non-smoothed) beam input, as a result, the final gain 
of it is smaller than the one computed with the initially 
smoothed beam. Fig. [8] shows the analytically-computed 
gain curve of the smoothed beam; Fig.[3]and Fig.fTUlshow 
the uncorrelated (slice) energy spread of the smoothed 
and the non-smoothed beam before BC1 and BC2 by 
ELEGANT, respectively. All the parameters in the ana- 
lytical calculation are obtained from ELEGANT, as what 
we do to the green curve in Fig. [7] The term "smoothed 
beam" means the whole 6 dimensional beam distribution 
is smoothed by removing the noise fluctuation, which is 
done by the sdds command smoothDist6s [12j . With this 
command, the noise fluctuation of the beam current pro- 
file in both the longitudinal and the transverse directions, 
and that of the uncorrelated (slice) energy spread are all 
suppressed to almost zero. We can see in Fig. [5] that the 
peak of the curve is much higher than that obtained from 
the noisy beam due to the smaller slice energy spread, 
and the wavelength where the peak resides is also much 
shorter, i.e., f - 20/im. However, although the total gain 
is small for the non-smoothed beam, the microbunching 
instability still can be a problem because of the signifi- 
cant initial noise. As the demonstration, Fig. [TT] shows 
that the final current fluctuation of the non-smoothed 
beam is comparable to that of the smoothed one, which 
is also consistent with the analytical results (Fig. [7] and 
Fig. [5J in terms of the wavelength where the peak re- 
sides. In summary, careful studies are needed to decide 
how much noise should be included in the computation. 
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FIG. 8: (Color) The gain curve calculated with the smoothed 
beam current. 
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FIG. 9: (Color) The slice energy spreads of the smoothed 
(red) and the non-smoothed beam (black) before BC1. 
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FIG. 10: (Color) The slice energy spreads of the smoothed 
(red) and the non-smoothed beam (black) before BC2. 
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FIG. 11: (Color) The current spectra of the smoothed (red) 
and non-smoothed beam (black) at the exit of the SXFEL 
linac, the wavelength is reduced by the factor of 10 from the 
initial value because of compression. 

The wavelength and the magnitude of the peak ap- 
pear to be slightly different in Fig. [7] for the two different 
methods. Possible reasons are: 

I. Beam parameters. As described, the beam param- 
eters used in the analytic formulae are obtained by EL- 
EGANT. Because ELEGANT and IMPACT-Z employ differ- 
ent algorithms to compute rf field, space charge, etc., the 
beam parameters calculated by the two codes could be 
slightly different. On the other hand, since the gain curve 
is very sensitive to the beam parameters such as the lo- 
cal (slice) energy spread, the transverse size of the beam, 
etc., therefore small variations of them can change the 
gain of the instability appreciably. As a demonstration, 
we see in Fig. [TJ] that, as the local energy spread rises, 
the peak gain of the instability falls rapidly. We may 
therefore conclude that in reality the gain can deviate 
from what we expect because of the possible uncertainty 
of the beam parameters. 

II. Numerical smoothing. Usually when a particle-in- 
cell (PIC) code is executed, smoothing is performed au- 
tomatically for the particle density distribution in the 
computation of the space-charge fields because of the 
grid-based algorithm. The modulation of particle den- 
sity could be suppressed excessively depending upon the 
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FIG. 12: The peak gain of microbunching instability as a 
function of the local (slice) energy spread before the first com- 
pressor (BC1) in the SXFEL linac. 



choice of the mesh size. For example, if the mesh size is 
too large, the final gain computed by IMPACT-Z can be 
smaller than the realistic and the modulation wavelength 
will be longer. Moreover, for a formula-based simulation 
code like ELEGANT, smoothing is carried out by low-pass 
filters Q when wakefields are computed. Therefore the 
same problem also exists in ELEGANT. Example relevant 
to noise smoothing will be shown in the next section. 

III. Plasma effect. If density fluctuation exists, charged 
particles will move between the higher-density region 
and the lower-density region, which is known as space- 
charge (density) oscillation. It usually happens in low- 
energy high-intensity beams. The mechanism is similar 
to plasma oscillation [13J, [l4[ • Some works have already 
been performed in this direction [15| . Preliminary in- 
vestigation by solving the Vlasov and Poisson equations 
in the six-dimensional phase space shows that the LSC 
impedance (Eq. (|2j)) will be modified when the beam ve- 
locity information and the space-charge oscillation fre- 
quency are included [16j . In our case, the longitudi- 
nal plas ma oscillation wavelength calculated based on 
k p = ^2I e /( r yllA&x) P3 i s about 120 m, where k p is 
the wave number of the oscillation, I e is the electron 
peak current, 70 is the electron energy, a x is the hori- 
zontal beam size and I a ~ 17 kA is the Alfven current. 
Because the plasma wavelength is much larger than the 
length of the first acceleration section (LI) (rs 6 m) be- 
fore BC1, we thus do not need to pay too much attention 
on the plasma effect at this point. However, the plasma 
effect in the injector must be investigated because of the 
low electron energy. 



III. CRITICAL ISSUES IN THE STUDY 



A. Higher order harmonics of beam current 

Equation ([1]) only takes into account the gain of the 
fundamental mode of beam current. In other words, 
the gain curve obtained from it merely describes the be- 
havior of the fundamental mode and nothing else, e.g., 
Fig. [7] and Fig. [8] in this paper. The linear theory does 
not include the higher order harmonics, which may ex- 
ist all the time simultaneously. However, in our simula- 
tions, higher-order harmonics of beam current are also 
observed. We therefore believe that, in certain cases, 
higher harmonics should be included in the computation 
of the gain. Preliminary studies tell us that the gain 
of the higher harmonics depends not only on the peak 
current of the beam, but also on the modulation depth. 
Thus the inclusion of the higher-harmonics becomes a bit 
more complicated than the fundamental mode. As an ex- 
ample, Fig. 1131 shows the excitation of higher harmonics 
simulated by ELEGANT after the first bunch compressor 
(BC1) with initial modulation wavelength of 100 /im (cor- 
responding to 20 /im in Fig. I X3[) which is in the vicinity 
of the peak (Fig. [7]) . The initial modulation depth is 5% 
and The Fast Fourier Transform (FFT) spectrum of the 
initial beam current is illustrated in Fig. 1141 We see in 
Fig.[l3]that, besides the peak gain of the fundamental re- 
siding around 20 /im, there are also the first, second, and 
third harmonics excited at, respectively, 10, 6.7, and 5 
/im. Similar excitations are observed as well at other ini- 
tial modulation wavelengths such as 80 /im, 60 /im, and 
35 /ium. Thus we need to keep in mind that higher-order 
of the initial modulation can also be strongly excited and 
the linear growth theory should be extended. 
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FIG. 13: The Fast Fourier Transform (FFT) spectrum of the 
beam current with initial modulation wavelength of 100 fiva 
at the exit of BC1, where 20 /mi in this figure corresponds to 
the initial 100 /im because of compression ratio of 5. 



In our study, some new issues are found to have pro- 
found effects on the final results of the microbunching 
instability and they need to be included in the compu- 
tation. These issues include higher harmonics of beam 
current, transverse-longitudinal coupling in a bunch com- 
pressor, and the method of noise smoothing. In this sec- 
tion, we will discuss these issues one by one in detail. 



B. Transverse-longitudinal coupling 

One important topic in microbunching-instability 
study is noise. Like the self- amplified spontaneous emis- 
sion (SASE) process, micro-bunches can be generated 
from shot noise. If the beam is ideally smoothed in 




20 40 60 80 100120140 
wovelength (micron) 

FIG. 14: The Fast Fourier Transform (FFT) spectrum of the 
beam current with initial modulation wavelength of 100 /im 
at the exit of the injector. 



all directions, i.e., without any noise, fluctuation, mod- 
ulation, etc., microbunching instability can hardly be 
excited. To reduce microbunching, there are ways to 
smooth the longitudinal current profile from the cath- 
ode by improving the temporal stability of driving laser. 
However, our investigation indicates that the density 
noise in the transverse direction can also introduce rough- 
ness to the longitudinal beam density distribution as a 
result of transverse-longitudinal coupling, for example, 
the -R51 elements in the transfer matrices of the dipoles 
in a bunch compressor. 
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FIG. 15: (Color) The current (black) and energy (red) spectra 
of a beam smoothed in all directions at the exit of BC1 with 
the laser heater turned off. 

Assuming a beam smoothed in all directions, and an- 
other beam with ~ 5% noise level in transverse but 
smoothed longitudinally, Figs. [15] and [16] show the cur- 
rent and energy spectra of the beams at the exit of the 
first chicane (BC1) without laser heating. In the figures, 
one can see that the levels of density and energy fluc- 
tuation (modulation) of the 3-D smoothed beam are un- 
doubtedly smaller than those when the beam is smoothed 
only longitudinally. 

Nowadays, a laser heater has been a common device to 
suppress the microbunching instability 4] in the linac of a 
FEL facility. As an example, Fig. [T71 shows the structure 
of the SXFEL laser heater including an injection laser, 
an undulator, and a chicane-type bunch compressor. Be- 
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FIG. 16: (Color) The current (black) and energy (red) spectra 
of a beam smoothed only in longitudinal direction at the exit 
of BC1 with the laser heater turned off. 



cause of the transverse- longitudinal coupling, some ex- 
tra noises are introduced into the longitudinal direction 
and the amplitudes of the modulations are also ampli- 
fied when beam passing through the laser heater. Fig- 
ures [18] and [19] show the current and energy spectra of a 
beam longitudinally smoothed but transversely noisy at 
the location right after the designed location of the laser 
heater with and without the laser heater included in the 
simulations. We can clearly see that the laser heater 
introduces extra peaks at various wavelengths and also 
increases the amplitudes of the current and energy mod- 
ulation, which will compromise the beam quality there- 
after. For this reason, the -R51 element of the magnetic 
chicane should be revisited with more care. The param- 
eters of the SXFEL laser heater used in the simulation 
are: length 0.55 m, periods 10, undulator magnetic peak 
field 0.31 T, laser peak power 0.80 MW, laser spot size 
at waist 0.30 mm. 
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FIG. 17: The layout of the SXFEL laser heater at 130 MeV. 



C. Noise smoothing 

When computing wakefields, noise smoothing is always 
performed by the simulation code to obtain the beam dis- 
tribution. For example, in ELEGANT, the Savitzky-Golay 
filter is used to smooth the longitudinal beam current 
profile (thus reducing the noise) in the course of comput- 
ing various kinds of wakefields. Because the method of 
smoothing has a strong influence on the final results, it 
must be investigated with care. 
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FIG. 18: (Color) The current (black) and energy (red) spectra 
of a beam longitudinally smoothed but transversely noisy at 
the first BPM right after the laser heater in the SXFEL linac, 
with the laser heater turned on in the simulation. 
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FIG. 20: (Color) The local (slice) energy spread at the exit 
of the SXFEL linac with the laser heater turned on. The S-G 
filter half- width is set at 1 (black) and 10 (red). 
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FIG. 19: (Color) The current (black) and energy (red) spectra 
of a beam longitudinally smoothed but transversely noisy at 
the first BPM right after the location of the laser heater in 
the SXFEL linac, with the laser heater substituted by a drift 
in the simulation. 



The Savitzky-Golay (S-G) smoothing filter is a filter 
that performs a local polynomial regression (of degree k) 
on a series of values (of at least k + 1 points which are 
treated as being equally spaced in the series) to determine 
the smoothed value for each point [l8[. This filtering al- 
gorithm was first described in 1964 by Abraham Savitzky 
and Marcel J. E. Golay j5|. They proposed a method of 
data smoothing based on local least-square polynomial 
approximation and showed that fitting a polynomial to 
a set of input samples and then evaluating the resulting 
polynomial at a single point within the approximation 
interval is equivalent to the discrete convolution with a 
fixed impulse response. The low-pass filters obtained by 
this method are widely known as Savitzky-Golay filters. 
In general, there are two key parameters in this method, 
one is the order, the other is the half-width. The former 
represents the highest order in the polynomial fit, while 
the latter corresponds to the number data points needed 
to be sampled in the smoothing. 

In our study, we vary the half width of the S-G filter 
in the ELEGANT input, and discover that the difference 
between the output slice energy spreads with different 
S-G half-widths is significant (Fig. EH)) . This is under- 
standable because stronger suppression of noise results 



in stronger reduction of local energy spread (excessive 
smoothing). Therefore here comes the problem: how 
to differentiate the real noise from the numerical noise 
in simulation? We suggest to use the S-G filter with 
caution (e.g., using filter half- width < 5) and to use as 
many macro-particles as possible within the capability of 
hardware. In our ELEGANT simulations here, we fix the 
S-G filter half- width at 1 and employ 1 million macro- 
particles. 



IV. CONCLUSIONS 

The study of microbunching instability in the SXFEL 
linac is performed in detail as a classic example. The 
potential problems and some new effects are uncovered 
including the initial beam distribution, the higher or- 
der modes in the noise amplification, the transverse- 
longitudinal coupling and the plasma effect in the compu- 
tation of LSC. More investigations are needed to obtain 
deeper insights of those problems. 

Computations using both analytic expressions and 
numerical simulations show that the gain of the mi- 
crobunching instability indicates large discrepancy be- 
tween the noisy and the smoothed beam input. Since 
in the real case, there is always noise introduced by the 
shot noise, the laser power jitters, etc., one should con- 
sider how much noise to be included in the initial input. 
On the other hand, our work shows that in terms of the 
final density /energy fluctuation (or the bunching factor), 
the difference between the noisy and the smoothed input 
is not large. 

A particle-in-cell (PIC) code appears to be a better 
tool for the study of microbunching instability because 
of its intrinsic advantage — the electromagnetic fields are 
computed by solving fundamental physical equations, 
i.e., the Poisson solver. As a result, a PIC code can 
be better in revealing the detailed physical process of 
the instability than a code based on analytic formulae. 
However, like all the other codes, the noise-smoothing 
method must be investigated to see whether it will intro- 
duce extra uncertainty. The challenging problem is how 



to separate and discard the numerical noise from the real 
one. Currently, we tend to preserve the information in 
the beam as much as possible by employing the least 
smoothing. 

The higher-order harmonics of beam current that re- 
vises the linear theory in the estimate of the gain should 
also be included in the development of the analytic for- 
mulae, especially in the vicinity of the peak gain. The 
transverse-longitudinal coupling in a magnetic chicane 
will transport shot noise from the transverse dimension 
to the longitudinal and needs to be considered in the de- 
sign of the laser heater. In this sense, we believe that the 
microbunching instability in an electron linac cannot be 
completely avoided unless the beam is ideally smooth in 
the whole 6 dimensional phase space. 

Because both the analytic and numerical methods ex- 
hibit limitations in the estimate of the microbunching in- 
stability, systematic experimental measurements are de- 



sired to provide solid and full understanding of the physi- 
cal essence. Work is ongoing at the Shanghai Institute of 
Applied Physics (SINAP) to prepare these experiments. 
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